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K*" . Abstract 
00 

^ ! We propose and examine a new set of benchmark supersymmetric scenarios, some of 

O ■ which have non-universal Higgs scalar masses (NUHM) and others have gravitino dark matter 

O ' (GDM). The scalar masses in these models are either considerably larger or smaller than the 

i-C \ narrow range allowed for the same gaugino mass mi/2 in the constrained MSSM (CMSSM) 

JL ! with universal scalar masses mn and neutralino dark matter. The NUHM and GDM models 
Ph. 

^ . with larger mo may have large branching ratios for Higgs and/or Z production in the cascade 
decays of heavier sparticles, whose detection we discuss. The phenomenology of the GDM 
\ models depends on the nature of the next-to-lightest supersymmetric particle (NLSP), which 
^ \ has a lifetime exceeding 10^ seconds in the proposed benchmark scenarios. In one GDM 
scenario the NLSP is the lightest neutralino Xi ^ind the supersymmetric collider signatures 
are similar to those in previous CMSSM benchmarks, but with a distinctive spectrum. In the 
other GDM scenarios based on minimal supergravity (mSUGRA), the NLSP is the lighter 
stau slepton fi, with a lifetime between ~ 10"^ and 3 x 10^ seconds. Every supersymmetric 
cascade would end in a fi, which would have a distinctive time-of-flight signature. Slow- 
moving fi's might be trapped in a collider detector or outside it, and the preferred detection 
strategy would depend on the fi lifetime. We discuss the extent to which these mSUGRA 
GDM scenarios could be distinguished from gauge-mediated models. 
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1 Introduction 



Among the most promising prospects for new physics that might be discovered at the LHC is 
supersymmetry. Its appearance at the TeV energy scale is motivated by the stabihty of the 
mass hierarchy [1], the hghtness of the Higgs boson as inferred from precision electroweak 
measurements [2], the possibihty of unifying the gauge interactions [3], and (assuming that 
i?-parity is conserved) the existence of a natural candidate for cold dark matter (CDM) [4]. 
However, despite these phenomenological hints and the intrinsic beauty of supersymmetry, 
there is no direct evidence for its existence. Evidence for new physics at the TeV scale 
might be provided by the anomalous magnetic moment of the muon, — 2, if it is con- 
firmed that measurements are significantly different from the Standard Model prediction [5]. 
Supersymmetry could be one of the possible interpretations, but not the only one. 

In advance of experiments at the LHC and other accelerators intended to explore the 
TeV energy range in more detail, it is valuable to understand the variety of possible signa- 
tures of supersymmetry, taking into account the constraints imposed by present accelerator 
experiments as well as astrophysics and cosmology. Experimental signatures of supersym- 
metry depend quite sensitively on the possible sparticle masses, which in turn depend on 
the unknown mechanism of supersymmetry breaking. There have already been a number 
of surveys of the supersymmetric parameter space [6, 7, 8, 9], some focusing on specific 
benchmark points that exemplify distinct possibilities, and others tracking the phenomeno- 
logical variations along lines in this space that are consistent with the cosmological and other 
constraints. 

Most such studies have focused on the CMSSM, in which the supersymmetry-breaking 
gaugino, scalar and trilinear parameters mi/2,''^o ^-nd ^Oj respectively, are each assumed to 
be universal at some input GUT scale ~ 10^^ GeV, and the gravitino is assumed not to be the 
lightest supersymmetric particle (LSP). However, this is not the only possibility: universality 
is not strongly favoured within our current understanding of possible underlying theoretical 
frameworks such as string theory, and the gravitino might be the LSP. Therefore, in this paper 
we investigate benchmark scenarios that sample models with patterns of supersymmetry 
breaking different from that in the CMSSM. 

The most dubious of the CMSSM universality assumptions may be that for the soft 
supersymmetry-breaking scalar masses itlq. Universality between sparticlcs in the same 
gauge multiplet is inevitable, and universality between sparticlcs in different generations that 
share the same quantum numbers is motivated by constraints from flavour-changing neutral 
interactions [10]. Moreover, sparticlcs with different quantum numbers may originate from 
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common GUT multiplets, in which case their soft supersymmetry-breaking scalar masses 
should also be universal at the GUT scale, apart from possible non-universal GUT D terms. 
However, none of these arguments give any reason why the soft supersymmetry-breaking 
scalar masses of the electroweak Higgs multiplets should be universal, and this may be the 
Achilles heel of the CMSSM. Accordingly, some of the non-CMSSM scenarios we study in 
this paper have non-universal Higgs masses (NUHM) [11, 12]. 

We assume in this paper that i?-parity is conserved, so that the LSP is stable and a 
candidate for the CDM required by astrophysics and cosmology [4]. However, even if one 
accepts the CMSSM universality assumptions, it is not evident that the LSP is necessarily the 
hghtest neutrahno x, as usually assumed in CMSSM analyses. Another plausible candidate 
is the gravitino [13] - [16], the supersymmetric partner of the graviton, in which case options 
for the next-to-lightest supersymmetric particle (NLSP) include the hghtest neutrahno and 
the lighter stau slepton fi. The gravitino mass m3/2 is poorly constrained by accelerator 
experiments, and the astrophysical and cosmological constraints on the gravitino differ from 
those on the neutrahno. One approach to gravitino dark matter (GDM) models is to retain 
the standard CMSSM universality assumptions and simply assume that the gravitino is 
lighter than the lightest neutrahno x- One of the GDM benchmarks we explore is of this type, 
with a neutrahno NLSP. However, allowing for flexibility in the gravitino mass introduces 
an additional arbitrary parameter, and it is difficult to scan and characterize the larger- 
dimensional parameter space opened up in this way. 

As an alternative scenario with a lower-dimensional parameter space, we consider a 

sample of GDM models in which the gravitino mass is fixed to equal the universal soft 

supersymmetry-breaking masses of observable scalar sparticles: m2/2 — ttio at the GUT 

scale. This assumption is motivated by minimal supergravity (mSUGRA) models of su- 

persymmetry breaking, which also impose a relation between the soft trilinear and bilinear 

supersymmetry-breaking parameters: Aq — Bq + mo [17]. This latter assumption may be 

used to fix the ratio of MSSM Higgs vacuum expectation values tan /3, further reducing the 

GDM parameter space [18]. For definiteness, we use the value Aq — (3 — -\/3)mo found in the 

specific Polonyi model of supersymmetry breaking in a hidden sector [19]. Even with this 

extra assumption, there is still a two-dimensional region of the {mi/2,mo) parameter space 

allowed by the accelerator, astrophysical and cosmological constraints. Below we propose 

GDM benchmark scenarios that explore some of phenomenological possibilities in this case, 

in which the NLSP is the fi ^ 

^One could also discuss GDM models with non- universal Higgs masses, but this would increase the 
dimensionality of the parameter space still further. 
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The structure of this paper is as follows. In Section 2 we give an overview of the CMSSM, 
NUHM and GDM parameter spaces, and in Section 3 we discuss in more detail the proposed 
new benchmark scenarios, specifying the parameter choices and presenting their spectra as 
computed in the SSARD [20] and ISASUGRA [21] programs. We also present calculations of the 
LSP relic density [22, 23], the 6 — > 57 decay branching ratio [24] and the supersymmetric 
contribution to g^^u — 2 [5] in each of these scenarios, as well as the value of found in a 
global fit to laboratory observables including my^ and sin^^^K [25]. Next, in Section 4 we 
present and discuss several of the most important sparticle decay branching ratios in the 
various new benchmark scenarios. Subsequently, in Section 5 we discuss the detectability of 
these benchmark scenarios at the LHC [26, 27], ILC [28] and CLIC [29]. Section 6 contains 
a dedicated discussion of possible measurements of the metastable fi NLSP mass and its 
decays in the various mSUGRA-motivated scenarios. Finally, Section 7 discusses our main 
conclusions and possible directions for future work. 

2 Overview of the CMSSM, NUHM and GDM Param- 
eter Spaces 

Before discussing specific benchmark scenarios, we first give an overview of the CMSSM, 
NUHM and GDM parameter spaces, which helps to motivate the specific parameter choices 
we make. In both the NUHM and GDM models, we implement the LEP limits on m^± [30] 
and rrih [31] and the b ^ constraint [24] in the same way as already documented for the 
CMSSM [7, 9, 22]. We restrict our attention to models with > 0, as favoured by g^^u — 2 [5]. 
We do not impose a numerical range on the supersymmetric contribution to this quantity, 
but we do refer later to a likelihood analysis that incorporates it as well as niw, sin^ 9w and 
6 — > 57 in the global function [25] . 

We first consider the CMSSM, assuming that the gravitino is not the LSP. We recall that 
the cosmological GDM density constraint: 

0.094 < ^cDMh^ < 0.129 (1) 

provided by WMAP [32] and earlier data restricts the CMSSM parameter space to narrow 
strips in the (mi/2, mo) planes for specific choices of tan (3 and Aq, if one assumes that all the 
GDM is provided by the neutralino LSP. The WMAP strips foliate the (mi/2,mo) plane in 
the CMSSM, as seen in panel (a) of Fig. 1, where we display updated WMAP strips for fj, > 
0,Ao — and tan (3 — 10, 20, 35, 50 and 55, as calculated assuming rrit — 172.7 GeV, the 
current central value [33]. The lighter (darker) parts of these strips are (in) compatible with 
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Updated CMSSM benchmark scenarios 



Model 


A" 


B" 


C" 


G" 


H" 


r 


J" 


L" 


M" 


mii2 


600 


250 


400 


375 


910 


350 


750 


450 


1075 










(935) 








(1840) 


mo 
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(60) 


(85) 


(115) 


(245) 


(175) 
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tan/3 
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20 


20 


35 


35 


50 


55 
(50) 



Table 1: We list updated GUT-scale input parameters for CMSSM benchmark scenarios with 
II > in the coannihilation and rapid- annihilation funnel regions. In scenarios where the 

values o/mi/2,mo and/or tan f3 have been changed so as to keep Vt^h"^ within the WMAP 
range, as calculated with the latest version o/SSARD, the previous input values are listed below 
in brackets. These points are valid for nif = 172.7 GeV. 

5(^ — 2 at the 2-a level, if one uses e'^e" data to calculate the Standard Model contribution [5]: 

Sa^ = (25.2 ± 9.2) x 10"^°. (2) 

Also shown are updated versions A", B", ... of the previously-proposed CMSSM benchmark 
scenarios [9] that lie on these WMAP strips. The update from primed to double-primed 
points is largely due to the change in our adopted top quark mass (from 175 to 172.7 GeV [33]) 
and improvements to SSARD which now include full two-loop running of the RGEs. We recall 
that most of the regions below the WMAP strips are forbidden, because there the LSP would 
be charged, namely the fi, which would be stable in this case. 

As noted above, the latest version of SSARD has been used to calculate these WMAP 
strips. It has an improved treatment of higher-order corrections, which are significant at 
large tan (3 and particularly in the very sensitive rapid-annihilation funnel region, which is 
now barely visible for tan/3 = 50 and rut — 172.7 GeV, but present for tan/3 = 55, as seen 
in Fig. la. We have adjusted minimally the parameters of previous GMSSM benchmark 
points with /x > 0, so as to leave unchanged the previous collider phenomenology and keep 
the relic density fi^/i^ within the WMAP range. The corresponding changes in (mainly) mo 
and (occasionally) mi/2 and tan/3 are tabulated in Table 1 where we indicate the values of 
mi/2 and mo for the previous (single-primed) points [9]. Note that we have not updated 
the previous 'focus-point' benchmark scenarios (these are extremely sensitive to m^, whose 
measurements at the Tevatron collider are still evolving significantly), nor the points with 
II < that are disfavoured hy g^ — 2. 

For any given value of mi/2 and mo, the NUHM has two additional parameters, reflecting 
the degrees of non- universality 6u,d of the masses of the two MSSM Higgs doublets [11]: 
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Figure 1: The allowed regions in the (mi/2,'n^o) planes for (a) the CMSSM with a neutralino 
LSP [light/dark (blue) shaded strips], (h) the NUHM with a neutralino LSP for /i = 375 GeV 
and tua = 265 GeV [light (blue ) shaded strips], (c) the GMSSM with a gravitino LSP weighing 
10 GeV, and (d) mSUGRA with a gravitino LSP [light (yellow) shaded regions labelled r < 1], 
for A = Amo : A = 3 — Vs. In each case, the proposed benchmark scenarios are indicated by 
(red) crosses. Panel (a) displays WMAP strips for tan P = 5,10,20,35,50,55, panels (b, c) 
have tan/5 = 10, and in panel (d) tan/3 is fixed by the electroweak vacuum conditions with 
the values indicated along the (solid black) contours. 
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— (1 + ^u,d)nT'o- One representative example of an NUHM plane is shown in panel (b) 
of Fig. 1, where we have assumed for definiteness that /i — 375 GeV and ttia — 265 GeV, 
using rrit — 178 GeV ^. The pale (turquoise) region allowed by cosmology includes a small 
'bulk' region at low mi/2 and mo that extends into a short coannihilation strip. There is 
then a near-vertical rapid-annihilation funnel whose location is determined by the choice of 
rriA, which is flanked by very narrow allowed WMAP strips at both larger and smaller values 
of mi/2- There is then a continuation of the coannihilation strip and finally a third, broader 
vertical band where the relic density falls within the range allowed by WMAP. The latter 
band is a result of the fact that the LSP is becoming more Higgsino-hke as mi/2 > The 
dark (brick) shaded region at large mi/2 and small mo is forbidden because here the LSP 
would the (stable) fi and the medium (green) shaded region at large mi/2 and mo is excluded 
by 6 — > 57. Only regions between the two near- vertical black dot-dashed lines have effective 
potentials that are stable up to the GUT scale, and are hence permissible theoretically. 
The near-vertical black dashed and red dot-dashed lines represent the LEP constraints on 
m-x±,mh respectively ^. Thus, only the two narrow WMAP strips above the dark (brick) 
shaded region between the near- vertical red and black dot-dashed lines are consistent with 
all the constraints. 

We see that one of the options opened up by the NUHM is a range of values of mo that 
are considerably larger than the very narrow range of values allowed by the CDM constraint 
in the CMSSM coannihilation strip, for any fixed values of mi/2,tan/3 and Aq. In the next 
Section, we exploit this freedom to increase mo by proposing three larger-mo scenarios shown 
as red crosses in the (mi/2,mo) plane. One of these has the same values of /IjUIa as those 
chosen in panel (b) of Fig. 1, whereas the other two (indicated in brackets) have different 
values of /i, m^- As discussed in more detail below, the specific values of /i, niA, mi/2 and mo 
were chosen so as to offer various different sparticle cascade decay signatures including decays 
of the second-hghtest neutrahno X2 — hx, Zx that are not favoured in CMSSM scenarios 
with a neutralino LSP, as was discussed in [9] . 

Similar large values of mo are also attainable in the CMSSM, if the LSP is the gravitino. 
A representative sample (mi/2,mo) plane in such a GDM scenario is shown in panel (c) of 
Fig. 1, where we made the particular choice m3/2 = 10 GeV ^. The light (yellow) shaded 
region labelled by r < 1 is that allowed not only by accelerator constraints (the LEP m^± 

^At this value of tan (3 = 10, the NUHM plane for rrit = 172.7 GcV would be virtually identical. 

•^The elliptical solid black lines bound the preferred range of — 2, if the Standard Model contribution 
is calculated using e+e~ data alone. 

'^As in the case of panel b), we use here rrit = 178 GeV, but the plane for rrit = 172.7 GeV would again 
be virtually identical. 
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and nih limits are shown as near-vertical dashed black and dot-dashed red lines, respectively) 
but also by astrophysical and cosmological constraints [34] . Only below the diagonal dashed 
purple line can one satisfy the CDM constraint on the relic density of gravitinos produced 
in decays of the NLSP. However, the most stringent cosmological constraints in this GDM 
scenario come from comparing the baryon-to-entropy ratio inferred from measurements of 
the Cosmic Microwave Background (CMB) with that inferred from the measured Big-Bang 
Nucleosynthesis (BBN) calculations and hght-element abundances, which might have been 
altered by NLSP decay products [34, 15] ^. The light elements whose abundances we include 
in this analysis are ^He, ^He, Deuterium, ^Li and ^Li. Also shown in panel (c) of Fig. 1 is the 
dotted line where — mf^, above which the NLSP is the lightest neutralino x, and below 
which the NLSP is the lighter stau slepton fi, which decays predominantly into r gravitino. 
For comparison, the region with pale (turquoise) shading is the strip in the (mi/2, mo) plane 
that would have been allowed if the lightest neutrahno were the LSP. 

The GDM benchmark point 5 shown in panel (c) of Fig. 1 again has much larger mo and 
important X2 hx, Zx decays. In this case, since the NLSP is the x which has a lifetime 
~ 1.8 X 10^ s, it has no collider signature apart from missing transverse energy. This scenario 
therefore looks quahtatively similar to the CMSSM scenarios discussed earher [7, 9], apart 
from its different and larger value of mo. 

In probing the possibilities for GDM with values of mo below the m^ = nif^ line in 
the (mi/2,mo) plane, we propose below to study more restricted mSUGRA scenarios with 
m3/2 = mo, Aq — {Bq + l)mo and the Polonyi choice ^o = (3 — -\/3)mo [19], as seen in panel 
(d) of Fig. 1, which was produced using m^ = 178 GeV ^. The dashed, dot-dashed and dotted 
lines and the pale (turquoise) strip have the same significances as in panel (c). Also shown 
as solid lines are contours of tan /3, as fixed by the electroweak vacuum conditions. In the 
low-mo region, the NLSP is the fi, and its hfetime within the light (yellow) shaded region 
varies between ~ 3 x 10^ and 10^ s. 

We propose to survey this wedge-shaped GDM region by studying the three indicated 

points e, C, 77 located at the vertex and along the top and bottom sides of the wedge. In all 

these benchmark scenarios, the metastable fi would be detectable as a charged particle with 

an anomalously long time of flight. As pointed out in [35, 36], one might hope to trap some 

of the slower-moving charged NLSPs and detect their decays. However, the strategies and 

prospects for detecting and trapping fi's would be rather different for NLSPs with hfetimes 

^Here and in panel (d), we restrict our attention to regions of GDM parameter space where the NLSP 
lifetime tnlsp > 10'' s, for which the only relevant decay products are photons and electrons. 

^Lowering mt to 172.7 GeV would mostly affect the tan/3 contours, typically increasing tan/3 by ~ 4. 
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Figure 2: Left panel: Illustration of relevant aspects of the sparticle spectrum ordering and 
important decay modes in the (mi/2,mo) plane of the CMSSM for taxi (3 = 10,/i > and 
= 0. Right panel: Decay modes of the second neutralino X2 ^ (plus signs) and xh 
(crosses) in various NUHM and GDM benchmark scenarios, compared with the corresponding 
branching ratios along the the WMAP line for tan /5 = 10 and fi > discussed in [9]. These 
plots were obtained with PYTHIA 6.215 [37] interfaced to ISASUGRA 7.69 [21]. 

measured in hours or months, as we discuss later. 

3 Proposed New Benchmark Scenarios 

We now describe in more detail the proposed new benchmark scenarios. As shown in 
Fig. 1(a), previous benchmarks were located along the WMAP lines in the CMSSM pa- 
rameter space where 0.094 < Vt^h"^ < 0.129, which are very narrow in mo for any fixed values 
of mi/2,tan/3 and A^. As seen in Fig. 2, the orderings of sparticle masses vary in impor- 
tant ways across the CMSSM (mi/2,''Tio) plane, with important implications for the allowed 
and dominant sparticle decay modes. However, within the CMSSM, points with larger mo 
generally have larger values of Q^h'^, and hence are cosmologically unacceptable ^. 

^We recall, however, the allowed focus-point region at very large mo, whose location is very sensitive to 
mt- The central value of m* has changed significantly as new Tevatron Run II measurements have been 
taken into account [33], and, awaiting its stabilization, we do not discuss the focus-point region further in 
this paper. For reference, we note here also that we assume m*^'^(mf,) — 4.25 GeV, and that Aq — 0, except 
for the mSUGRA GDM scenarios. 
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In particular, along the CMSSM WMAP strips and specifically in the previous benchmark 
scenarios, the branching ratios X2 hx, Zx are generally quite small. The decays X2 
dominate in this region, as illustrated in both panels of Fig. 2. This is rather exceptional, 
since, as seen in the left panel, there is a large region of the CMSSM parameter space where 
X2 xh decays dominate. These could play an important role in the discovery of the lightest 
Higgs boson h, as well as in reconstructing the sparticle cascade decays [26, 38]. There is 
also a band of low values of mi/2 and relatively large values of mo where X2 decays 
are important, which does not happen along the WMAP strips, as seen in the right panel 
of Fig. 2. Even more strikingly, in models with values of mo significantly smaller than the 
WMAP strip value for any given choices of mi/2, tan /3 and Aq, the LSP would be the fi, but 
such a charged LSP would be forbidden by astrophysics. However, scenarios with values of 
mo very different from those in the WMAP strips can be consistent with cosmology beyond 
the specific CMSSM framework considered in [7, 9] , as we now discuss. 

3.1 NUHM Benchmark Scenarios 

When one considers models with non-universal soft supersymmetry-breaking contributions 

to the Higgs masses (NUHM), the electroweak vacuum conditions no longer fix |;Li| and the 

pseudoscalar Higgs mass uia as they do in the CMSSM, so that both mo and the X2 ^ X 

mass difference can be increased in such a way that the decays X2 xh, Z are kinematically 

accessible but not competing decay modes into sleptons. We have chosen to exemplify the 

new phenomenology these decay modes open up by proposing consideration of three NUHM 

benchmark scenarios with relatively small sparticle masses that are beyond the reach of the 

Fermilab Tevatron collider, but would provide plenty of early physics opportunities for the 

LHC and substantial follow-up opportunities for the ILC ^. 

As seen in Table 2 and as a (red) cross in Fig. 1(b), NUHM benchmark a has tan /? = 10 

and mi/2 = 285 GeV, which is close to the value mi/2 = 250 GeV in the CMSSM benchmark 

B" (or SPSla [8]). However, as seen in panel (b) of Fig. 1, it has a larger value of mo = 

210 GeV, as compared to the CMSSM benchmark B" value mo = 65 GeV. This ensures that 

X2 decays are kinematically forbidden, as well as X2 X^, whereas X2 xZ is possible 

and prominent, thanks to the mass difference m^j — — 99 GeV. As also seen in Table 2, 

this benchmark has /i — 375 GeV and m^ = 265 GeV, as for the (mi/2, mo) plane shown in 

panel (b) of Fig. 1. The choice of ttia has the important consequence that rapid annihilation 

via the A pole occurs at lower mi/2, reducing the x LSP density to the allowed cosmological 

^Their values of mi/2,mo and tan/3 approximate those of benchmarks chosen for study by the CMS 
Collaboration [38]. 
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range, as seen from the location of point a on one of the thin WMAP strips in panel (b) of 
Fig. 1 ^. We also note the degrees of non-universality in the soft supersymmetry-breaking 
Higgs masses: 

m]j^ = -(333GeV)^ m]j^ = +(294GeV)^ (3) 

where Hu {H^) give masses to the w-type quarks (d-type quarks and charged leptons), respec- 
tively. These correspond to violations of Higgs universality that are 0{1) . Coincidentally, 
the value of mi/2 at the benchmark point a is close to the value of mi/2 at the tan/? = 10 
CMSSM point that minimized a fit to accelerator data including m\\/,iiiv? 6y^r,b — > 57 

and — 2 [25]. The value of is not very sensitive to mo, so its value at benchmark a is 
not very different from the minimum value in the CMSSM, as seen in Table 3. This Table 
also lists this point's values of O^/i^, 6 — > 57 and the supersymmetric contribution to f?^— 2. 

As also seen in Table 2 and as a bracketed (red) cross in Fig. 1(b), NUHM benchmark f3 
has tan /? = 10 but the somewhat higher value of mi/2 = 360 GeV and mo = 230 GcV. How- 
ever, it has /i = 500 GcV and niA = 325 GcV, for which choices the (mi/25'^0) plane would 
look somewhat different. As for point a, the choice of m_4 ensures a suitable relic density via 
rapid annihilation. The non-universal soft supersymmetry-breaking Higgs masses-squared 
are: 

= -(461GeV)2, m|^ = +(229GeV)2. (4) 

These again correspond to violations of Higgs universality that are 0{1). With the larger 
value of // as compared to point a, the X2 — X mass difference is increased to 133 GeV, 
thereby opening up the decay mode X2 xh {iTih = 117 GeV at this point). Benchmark j3 
lies somewhat further from the GMSSM point with minimum x^, but its overall value of 
is still acceptable, as seen in Table 3, along with its values of fl^h"^, b ^ s'-f and Qfj, — 2. 

As also seen in Table 2, the final NUHM benchmark 7 has tan/3 = 20 and, as seen as 
a bracketed (red) cross in Fig. 1(b), the somewhat lower value of mi/2 = 240 GeV, and the 
somewhat larger value of mo = 330 GeV. It also has fi = 325 GeV and m^ = 240 GeV. The 
latter choice ensures an acceptable rehc density, although with different values of fi and uia 
from those used in the rest of panel (b) of Fig. 1. The non- universal soft supersymmetry- 
breaking Higgs masses-squared are: 

ml^ = -(242GeV)2, m^^ = +(373GeV)2. (5) 
^There are similar effects for the other NUHM benchmarks introduced below. 

^^Despite the negative value of rnjf^, the contribution of \fi\^ to the effective Higgs masses ensures that 
there is no vacuum instability below the GUT scale, apart from that induced at the electroweak scale by the 
top quark. 
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Supersymmetric spectra in NUHM and GDM benchmark scenarios 
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Table 2: Proposed NUHM (a, /3, 7^ and GDM (6, e,(,v) benchmark points and mass spectra 
(in GeV), as calculated using SSARD [20] and FeynHiggs [39], using the one-loop corrected 
effective potential computed at the electroweak scale and one-loop corrections to the chargino 
and neutralino masses. We use here mb{mb)^^ — 4.25 GeV and rut — 178 GeV. 
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Relic density Q.^h?^ b ^ and — 2 in post-WMAP benchmark scencirios, also 

tnlsp in the GDM models 
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Table 3: Comparison ofD,LSph^ for the benchmark points in Table 2, as computed with the 
SSARD code [20], 5a^(^10~^j, the branching ratio for b s^, the NLSP lifetime for GDM 
scenarios and the for a global fit to precision observables [25]. 



These also correspond to violations of Higgs universality that are C(l), not orders of mag- 
nitude. With the choices of tan/5,mi/2 and fi, the X2 — X mass difference is decreased to 
83 GeV, closing off the X2 — ^ X^y h decay modes. Moreover, with this large value of mo, the 
sleptons are far from the kinematic range for X2 decays, which are therefore dominated by 
non-specific three-body decays mediated mainly by virtual Z exchange. As seen in Table 3, 
the overall x"^ is again quite small for this point, which has a value of mi/2 rather similar to 
CMSSM benchmark B" ^\ 

3.2 GDM Benchmark Scenarios 

As discussed in the Introduction, one may construct GDM scenarios with values of thq that 
are either much larger or much smaller than that in the CMSSM for the same values of 
tan/3,mi/2 and Aq. 

GDM benchmark S described in Table 2 has been chosen to exemphfy phenomenology at 

a significantly larger value of mo, specifically 500 GeV, as well as a moderately large value 

of mi/2 = 750 GeV. Thus, it samples the 'panhandle' reaching to large mo, as shown by the 

(red) cross in Fig. 1(c). The X2 — X niass difference increases to 302 GeV in this scenario, and 

all the sparticles are significantly heavier than in the previous NUHM benchmarks. We recall 

that the NLSP in this scenario is the hghtest neutralino x- However, the neutralino lifetime 

for decays into a gravitino is 1.8 x 10^ s, too long to be detectable in the neighbourhood of 

a plausible collider experiment. Because of its larger values of mi/2 and mo, in particular, 

the overall x^ for this model is somewhat further from the global minimum, but it cannot 

^^We note that the spectra for the NUHM benchmark points were computed using rrit = 178 GeV. The 
primary cflFcct of lowering m,i to 172.7 GeV would be to lower mh by 2-4 GeV: squark and slepton masses 
would only change by 0-2 %. 
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be excluded on this ground. 

Here and in the following GDM models, we obtain a contribution to the LSP rehc density 
by first calculating the NLSP abundance following thermalization, annihilation and freeze- 
out, and then reducing the density by a factor m3/2/mjvL5P to allow for the subsequent NLSP 
decays to the gravitino LSP. We see from Table 3 that this contribution to the relic gravitino 
density lies below the preferred WMAP range, and this feature is even more marked for the 
other GDM models discussed below. The total cold dark matter could be raised into the 
WMAP range either by another gravitino production mechanism, such as thermal production 
in the very early Universe [13], or if there is another component in the cold dark matter, 
such as an axion or superheavy relic. 

Benchmark scenario 5 is just one point in a three-dimensional space parametrized by 
(mi/2, mo, m3/2) for the specific choices Aq — and tan/3 = 10, which certainly includes 
regions that would be more accessible in the early days of the LHC or at the ILC. On the 
other hand, the high-mo panhandle extends to larger values of mo that are not shown. This 
region presumably also extends to larger values of mi/2 than those shaded in Fig. 1, but these 
would have lower NLSP lifetimes, in which cases hadronic NLSP decays would also have to be 
taken into account when evaluating the astrophysical and cosmological constraints [34] . We 
do not explore these options, since the challenging nature of benchmark scenario S already 
serves as an adequate counterweight to the 'easy' scenarios a, (5 and 7 with their relatively 
low values of mi/2 and mo. 

Since benchmark point 5 is CMSSM-like, changing rrit to 172.7 GeV would produce very 
httle change in the sfermion spectrum, the biggest effect being 1.5 % in mf. The Higgs mass 
would drop by 4 GeV, m^ would drop by 3 %, and would drop by 5 %. The latter 
two changes would also affect the heavy Higgses and heavy neutralinos and charginos by 
approximately 5 %. 

In order to reduce the dimensionality of the parameter space to be explored by the 
remaining GDM benchmark scenarios, we next assume an mSUGRA framework in which 
m3/2 = mo and Aq — Bq -\- mo- We further assume that ^Iq = (3 — -\/3)mo [19]. In this 
case, the value of tan/3 is fixed by the electroweak vacuum conditions and varies across 
the (mi/2,mo) plane, as seen in panel (d) of Fig. 1. In addition to a (pale blue) WMAP 
strip with x LSP at mi/2 < 1 TeV, we also see a (yellow) GDM region of parameter space 
allowed by the astrophysical and cosmological constraints. It takes the form of a wedge that 
broadens as mi/2 increases, throughout which the NLSP is the fi. We choose as our next 
GDM benchmark e the point shown as a (red) cross that is close to the vertex of this wedge, 
with mi/2 = 440 GeV and mo = 20 GeV. In this case, tan/3 = 15, the fi NLSP has a mass 
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of 150 GeV and a lifetime of 3.3 x 10^ s. At this point, the global is not far from the best 
fit in the CMSSM, as seen in Table 2. 

As seen by two more (red) crosses in panel (d) of Fig. 1(d), we complement this LHC- 
and ILC-friendly point with two points that are more challenging. A priori, values of mi/2 
considerably beyond 2 TeV would be possible in this wedge: they would be beyond the reach 
of either the LHC or the ILC, although the CLIC reach in fi pair production would extend 
beyond mi/2 = 4 TeV. We consider two points with mi/2 = 1 TeV, which are already quite 
challenging for the LHC. The upper edge of the wedge is defined by the astrophysical and 
cosmological constraint on fi decays, and corresponds to a fi lifetime ~ 3.3 x 10^ s. The 
lower edge of the wedge that we consider corresponds to a hfetime of 10^ s. Benchmark 
point C is close to the upper edge, with mo = 100 GeV and tan /3 ~ 21.5. Here the fi NLSP 
has a mass of 340 GeV and a lifetime of 2 x 10^ s. Benchmark point rj is close to the lower 
edge, with mo — 20 GeV and tan /3 = 23.7. Here the fi NLSP has a mass of 322 GeV and a 
lifetime of 6.8 x 10^ s. Both these points have rather larger values of than the best fit in 
the CMSSM, as also seen in Table 2, but these points cannot be excluded on these grounds. 

The changes in the spectra due to the shift in rut would be as follows: tan /3 would be 
increased by ~ 4, would be lowered by 4 % and rriA would be lowered by 5-7 %, with 
corresponding changes in the heavy Higgses, neutrahnos and charginos. The hght Higgs 
mass would be lowered by 3-4 GeV, and changes in the sfermion masses would typically be 
less than 1 %, with the exception of the lighter stau, whose mass would drop by about 7 %. 

3.3 Discussion of Spectra 

As in our previous papers on CMSSM benchmarks [7, 9], the parameters of the NUHM and 
GDM benchmarks were first specified using the code SSARD [20]. In order to facilitate the 
interfaces with standard simulation packages, the spectra calculated with SSARD were then 
matched using parameters of the ISASUGRA 7 . 69 code [21] to reproduce the main features 
of the SSARD spectra. The values of mo and mi/2 were adjusted to give the same masses 
for the lightest neutralino x and the lighter stau fi. Then the Higgs mass parameters ttihu 
and rriH^ were varied to reproduce the SSARD values of m^ and jj,. As these choices altered 
slightly the values of m^ and m^j , the procedure was then iterated. The final ISASUGRA 7 . 69 
parameters are hsted in Table 4. This procedure was not followed for the GDM points, as 
our results are less sensitive to the exact spectra, and here the SSARD inputs were used. Note 
the difference in the sign convention for Aq between the two codes. 

As already mentioned, the first three of the new benchmarks, a,/?, 7, are NUHM points 
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Supersymmetric spectra in NUHM and GDM benchmcirks 
calculated with ISASUGRA 7.69 
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Table 4: Proposed NUHM and GDM benchmark points and mass spectra (in GeV), as cal- 
culated using ISASUGRA 7.69 [21] and adapting the input parameters to give the best match 
to the SSARD [20] spectra shown in Table 2, as described in the text. 
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chosen to yield rather low-mass spectra, observable at an early stage of the LHC running, 
as might also point e. They also offer interesting physics opportunities for the ILC. These 
points complement the previous benchmark points B', C and I' of [7, 9], as they give rise 
to different search topologies. On the other hand, points S, ( and rj have heavier sparticles, 
and hence are much more challenging for both the LHC and the ILC. 

4 Sparticle Decay Branching Ratios 

One of the key particles appearing in sparticle decay chains is the second neutralino X2, 
whose branching ratios are quite model-dependent and have significant impact on sparticle 
detectability at future coUiders [6, 7, 8, 9, 26, 27]. In particular, X2 decays play crucial roles 
in reconstructing the masses of heavier sparticles such as squarks and gluinos via cascade 
decays. Moreover, X2 decays may offer new ways to discover or measure other new particles, 
such as the lightest MSSM Higgs boson h. Therefore, we now use ISASUGRA 7.69 to discuss 
the principal branching ratios of the X2 in the various NUHM and GDM benchmark scenarios 
introduced above, comparing them in particular with those in the CMSSM at different points 
along the WMAP hne for tan/3 = 10 and > 0, as discussed in [9]. 

We first recall the principal branching ratios of the X2 in the low-mass CMSSM bench- 
marks considered previously. In the case of point B', the dominant decay mode was X2 ^r^ 
(11 %) followed by £r whereas in the case of point C the dominant decay mode was 

X2 — (11 %) followed by £i — > x^- On the other hand, in the case of point I', which has 
a relatively large value of tan /3 = 35, the dominant decay was X2 tit (96 %), followed by 
fi XT- 

The new points a, /3, 7 provide qualitatively new signatures, as shown in Fig. 2. At the 
point a, the X2 mainly decays via X2 Zx (96 %), which is observable through the Z 
leptonic decay mode. At the point /3, the main decay signature is X2 hx (64 %), where 
the Higgs boson can be reconstructed from its decay to bb. In addition, there are smaller 
branching ratios for X2 — Zx (8 %) and X2 — ^ t^iT (23 %). On the other hand, point 7, 
which is just above the Higgs mass limit from LEP, has direct three-body leptonic decays 
X2 (-(-X (4 %) and X2 tt'X (3 %), and the other decays are mainly X2 QQX niediated 
by virtual Z exchange. 

At all three points, the chargino decays dominantly into Wx- At points a and /3, the 
gluino is heavier than any of the squarks and decays to qq. The qu decays directly to qx, 
whereas the qi leads to cascade decays such as qL — > qX2 (typically ~ 30 %) and qL — > q'x"^ 
(typically ~ 60 %). On the other hand, at point 7 the gluino is hghter than the squarks of 
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the first two generations, and its dominant decay is ^ — > bbi (81 %), followed by 61 — > bx2 
(26 %), 61 ^ tx^ (36 %) or 61 — > txf (26 %). The squarks of the first two generations decay 
similarly as at points a and (3 and with similar branching ratios. 

As already remarked, at point 5 the NLSP is the neutralino, which looks stable from the 
point of view of a collider detector, and gives rise to the usual missing-energy signature. As 
seen in Fig. 2, a further signature is provided by X2 — hx (91 %), with a smaller branching 
ratio for X2 Zx (5 %). At this point also, the gluino is heavier than any of the squarks, 
whose decays are similar to those at points a and /3. 

At the last three points eX,V with a fi NLSP, the gluino is heavier than any squark and 
decays to qq with some preference for iit (~ 20 %). The most important products of the 
subsequent squark cascade decays are displayed in Table 5, which we now explain. It is a 
general feature of these models that qji — > qx with branching ratios ~ 100 %. Then x — > fiT 
with large branching ratios of 92/75/69 % in models e, C,?7, respectively, with essentially 
all the other decays being x ^R^ followed by Ir — > Itrt. As a result, the dominant qR 
final states are qfiT, with somewhat smaller fractions of qUf\T. Analogously, in many cases 
Q.L — qX2-i with branching ratios of 32/33/33 % in models e,C,,r], respectively, the other 
decays mainly being qL — > q'x"^- Many of the subsequent X2 decays are also to fir, or else 
IrI followed again by Ir — > Itrt Thus, qi, decays via X2 populate the final states qfiT, 
qiifiT and qiii'i'fiT. On the other hand, qi q'x^ decays mainly populate final states 
containing neutrinos that would be more difficult to reconstruct, with the possible exception 
of some q'rWfi final states. 

The general conclusion is that LHC final states in these GDM models with fi NLSPs 
contain a pair of r leptons and quite possibly additional lepton pairs, as seen in Table 5. In 
the case of benchmark scenario e, — m^^ ~ 35 GeV and m-^ — vrif^ ~ 23 GeV, so the 
efficiency for picking up the additional cascade decay leptons may be reduced at the LHC 
but in benchmark scenarios C, V these mass differences exceed 100 GeV, and these cascade 
decay leptons should be readily detectable. 

^^There may also be some qL qf2T decays, followed by f2 fiZ/h decays. 

^•^Wc note, however, that it will not be necessary to trigger on these leptons, since the sparticle production 
events will generally contain many energetic hadronic jets. 
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Table 5: Final states for the benchmark points with Gravitino Dark Matter (GDM), as 
calculated with ISASUGRA 7.69 [21]. 

5 Observability at Different Accelerators 
5.1 Detect ability at the LHC 

We now provide rough estimates of the numbers of different species of supersymmetric parti- 
cles that may be detectable at the LHC in the various NUHM and GDM benchmark scenarios 
introduced above. We assume the production cross sections for squarks and gluinos at the 
LHC that are listed in Table 6. The physics objects shown in the figures below are obtained 
in the following way. 

• Jets arc reconstructed from particles generated by the PYTHIA [37] Monte Carlo, using 
an iterative cone algorithm with a cone size of 0.5 radians. In order to model a typical LHC 
detector acceptance, wc require each jet to have a pscudorapidity |?7| < 3.0 and a transverse 
energy Et > 20 GcV. These jets include hadronic tau decays. 

• The missing transverse energy is calculated from the transverse energies Et of the 
visible particles. 

• Charged leptons e, fi arc accepted if their transverse momenta Pt > 10 GcV and their 
pseudorapidities \ri\ < 2.4. Their momenta are smeared with a Gaussian error between 1 % 
and 10 %, depending on the momentum. 

• We assume a 50 % efficiency for identifying b jets, with mis-tagging rates of 15 % for 
charm jets and 5 % for light quarks and gluons. 
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• We assume a 50 % efficiency for identifying fiadronic r decays, witfi a 6 % mis-tagging 
rate for jets with Et < 30 GeV and a 1 % mis-tagging rate for jets with Ex > 30 GeV [26]. 
More complete and sohd results should be obtained from detailed experimental simulations. 
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Table 6: Cross sections in pb for models with non-universal Higgs masses (NUHM) or 
Gravitino Dark Matter (GDM) calculated with PROSPINO [40] at NLO and with masses 
from SSARD. The squark cross section is computed for the first 5 flavours. The quantity 
atot is the inclusive sum over all production mechanisms of the gluino or squark, e.g., 
<^tot{g) = 2 o'{gg)a{qg). Production of ii and ^2 ^re not included in the sums, but listed 
separately. 

We start with the benchmark points a to 5, adopting criteria similar to those used 
previously in discussions of CMSSM benchmark scenarios with a neutralino LSP [7, 9] . 

• Higgs bosons: We generally follow the ATLAS and CMS studies of the number of 
observable Higgs bosons as a function of tha and tan/3 [26, 27], bearing in mind that they 
have no significant 'exotic' decay modes into non-Standard Model particles. The hghtest 
neutral Higgs boson h is detectable at all four points, and the heavier neutral Higgs bosons 
H, A would be observable in scenarios a, /3 and 7. In contrast, the charged Higgs bosons 
i?^ would be observable only at point 7 where tan P — 20 since, according to previous 
studies [26, 27], if^ cannot be seen at the LHC when tan/3 = 10, for any studied values of 
the other MSSM parameters. 

• Gauginos: The lightest neutralino x is considered always to be observable via the 

cascade decays of observed supersymmetric particles We consider the X2 to be observable 

at the LHC if the product of its production cross section and the relevant decay branching 

ratio (xh, Z or i'^£~) is at least 0.01 pb, corresponding to 1000 events produced with 100 fb~^ 

of integrated luminosity. Thanks to the rather large production cross sections in Table 6, 

the X2 is observable at all four points in the cascade decays of squarks. On the other hand, 

^^We recall that, from the detector point of view, the x is effectively stable in benchmark 6, and hence 
has a missing transverse energy signature, as at points a, /3, 7. 



19 



as the lighter chargino decays with a branching ratio > 90 % into Wx, it will be difficult 
to detect in cascade decays, so the possibility we have considered is via direct production of 
X2X^, leading to tri-lepton final states. Previous studies have indicated that the x^ would 
not be observable in this mode for mi/2 > 170 GeV, as in all the benchmark scenarios 
a, P, 7, 5. On the other hand, we note that the associated X2X^ production cross sections at 
points a, 7 are ~ 2.5, 5.1 fb, respectively, so these points would be worth further study. 

One of the motivations for specifying the parameters of scenarios /3 and S was to consider 
models with large branching ratios for the cascade decays X2 X^, which have been studied 
previously by both ATLAS and CMS. Fig. 3 shows the cascade /i — > 66 signals expected in 
scenarios P and S. We select events with missing Et > 150 GeV and require the candidate b 
jets to be separated by an 'angle' AR = V Ajy^ + < 2. We have generated 5000 sample 
supersymmetric events in each scenario: these correspond to integrated luminosities of less 
than one fb~^ for scenario /3 and about 50 fb~^ for scenario S. The solid lines are the signals 
and the dashed lines are the supersymmetric backgrounds in the two scenarios from events 
not containing Higgs bosons the Standard Model backgrounds are much smaller. The h 
is visible in X2 decays in scenario /3 (we recall that the event sample corresponds to a very 
small integrated luminosity in this case) and the signal is even clearer at point 5. Previous, 
more detailed CMS and ATLAS studies of low-mass points similar to /3 have also concluded 
that the h should be observable in sparticle cascade decays, and a point very similar to /3 is 
being studied thoroughly for the CMS physics TDR [38]. Our first examination of point S 
is promising, but more detailed studies of such a high-mass point would also be useful. 

• Squarks: The spartners of the lighter quark fiavours d, s, c are considered to be 
observable if rriq < 2.5 TeV [41], so they could be observed in all four benchmark scenarios. 
However, their fiavours could not be distinguished at the LHC. We further assume that the 
stops and sbottoms f, b are identifiable only if they weigh below 1 TeV, unless the gluino 
weighs < 2.5 TeV and the stop or sbottom can be produced in its two-body decays. As 
in scenarios a, (3, 7 the stops and sbottoms are relatively light, we consider them to be 
observable. At point 6, the branching ratio for ^ — > fit ~ 40 %, whereas the decays into 
t2t, bib and 626 are each only (9(10) %. Accordingly, we consider only the ti to be detectable 
at point 5. A note of caution is in order: since the detection of f or 6 production is difficult 
to assess without simulation studies, these conclusions should be taken with care. 

We have considered the detectability of qr via their decays qr — > X2Q followed by X2 
Xh. Selecting the events in the h peaks in X2 decays shown in Fig. 3(a,b), and then combining 
the reconstructed h boson with the hadronic jet Ji that maximizes the product ph • pji , we 

These are mainly due to events containing b squarks. 
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Figure 3: The signals from the lightest Higgs boson h in X2 ^ cascade decays in scenarios 
(a) P and (h) 5, respectively, as found using PYTHIA [37] interfaced with ISASUGRA [21]. In 
each case, the solid line shows the sum of the super symmetric signal and background, and 
the dashed line is the supersymmetric background alone. The Standard Model background is 
much smaller. The qh invariant mass distributions in scenarios (c) (3 and (d) S, respectively, 
each exhibit an end-point corresponding to Qr qx2 cascade decays followed by X2 xh- 
The qqh invariant mass distributions in scenarios (e) [3 and (f) 6, respectively, each exhibit 
a feature corresponding to g ^ q^q and q^ qx2 cascade decays followed by X2 X^- 
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obtain the candidate qh invariant mass distributions shown in panels (c, d) of Fig. 3 for points 
/3, S, respectively. In principle, qr decays give end-points in these distributions, because of 
the two-body decay phase space, which are less distinctive than the corresponding dilepton 
edges in — > qt^l~x decays. Edge features are visible in both scenarios /3, 5, close to the 
expected positions at ~ 600, 1200 GeV, respectively 

• Gluinos: These are generally considered to be observable for masses below 2.5 TeV [41], 
and hence can be discovered at all four points. 

As seen in panels (e, f ) of Fig. 3, we have also considered specific features of the searches 
for gluinos at points /3, 5, respectively. In each case, we have selected the hadronic jet J2 that 
maximized the product pj^ • pjj and plotted the hJ\J2 mass distribution. Gluino decays 
should give distinctive 'edge' features at ~ 640, 1280 GeV, respectively, corresponding to the 
multi-body phase space. This feature is not very apparent for point (5 with the small sample 
generated here, but more apparent for point 5. A point similar to the former is also under 
study for the CMS Physics TDR [38] . A more detailed study of benchmark 5 with optimized 
cuts would also be desirable. 

We have also considered the decays ^ — > 66 at points q;,/3, 7 and 5. In each case, the 
product of the g production cross section and branching ratio appears high enough to enable 
rrig — mi to be measured with sufficient accuracy to verify that this point has a value of 
mo significantly different from that on the CMSSM WMAP hne for the same value of mi/2. 
However, more refined studies of gluino search strategies would clearly be useful. 

• Charged Sleptons: Since the mass differences — ~ 100 GeV at all these bench- 
marks, we consider i decays into leptons always to be observable, provided that the slepton 
mass is light enough and hence the production cross section is large enough. We note that 
all four points have negligible branching ratios for the decays of X2 to sleptons, which im- 
plies that cascades will not contribute to slepton observability and that sleptons can only be 
detected via their direct production. Following [7] , we consider the direct production rates 
to be large enough if < 350 GeV. According to this criterion, all the charged sleptons 
would be observable in scenarios a and /3 (though £l and f2 signals would be very marginal 
at the latter point), and cr, jiR, fi would be observable at point 7. However, we consider the 
observability of the fi^2 to be difficult to assess without a detailed study and, conservatively, 
we do not count them as observable in any scenario. The sleptons are all too heavy to be 
observed at point S. 

• Sneutrinos: We do not consider sneutrinos to be observable at the LHC. Although many 
of the sneutrino decays are into visible particles at points q;,/3, 7, the modes (A*^X^,i^X2) 

Similar distributions are being studied in more detail for the CMS Physics TDR [38]. 
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Figure 4: The missing Et spectra in scenarios e and C,: that for point rj is similar to the 
latter. These plots were obtained with PYTHIA [37] interfaced with ISASUGRA [21]. 

having branching ratios (46, 20), (37, 17), (55, 33) % respectively. Moreover, the associated 
iu production cross sections are quite large in scenarios q;,/5, 7: 220/110/80 fb, respectively. 
Nevertheless, no viable discovery strategy has yet been developed. 

We next discuss particle observabilities in the GDM scenarios e, (, rj, where fi is the 
NLSP. The branching ratios for final states resulting from the interesting cascade decays of 
squarks at the LHC are shown in Table 5. We first note that sparticle pair-production in 
these scenarios gives rise to substantial missing Et, as seen in Fig. 4 for scenarios e and ( 
(point rj is very similar to the latter). We assume that the metastable fi's are measured in 
the detector, as discussed below. The missing Et is traceable to the many neutrinos in the 
final states, e.g., from r decays and/or the many qi decays with other neutrinos. 

We display in Figs. 5, 6 and 7 some other characteristics of events in these benchmark 
mSUGRA GDM scenarios. Panels (a) of Figs. 6 and 7 show the jet multiplicity distributions 
for points (, rj. Two-body decays of the are responsible for the bimodal distibutions of the 
leading jet transverse energies in panels (b) of both Figs. 6 and 7. The peak at Et ~ 1 TeV 
is due to the two-body q^ decays, and the lower- peak is due to other sparticle decays. 
Panels (c) of Figs. 6 and 7 show barely visible features in the leading lepton Et distributions 
due to slepton cascade decays. There is no such feature in Fig. 5, where the cascade lepton 
energy is smaller. Nevertheless, we note that large fractions of the cascade-decay leptons have 
transverse momenta large enough to be detected with high efficiency, and could potentially 
be used as event triggers in addition to the high-ii^^ jets. 
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Figure 5: Kinematic distributions in events produced in the GDM benchmark scenario e: 
(a) hadron jet multiplicity, (b) the transverse energy Et of the most energetic jet, (c) the 
transverse energy pt of the most energetic lepton, and (d) the pseudorapidity distribution for 
fi production. These plots were obtained with PYTHIA [37] interfaced with ISASUGRA [21]. 
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Figure 6: Kinematic distributions in events produced in the GDM benchmark scenario (, as 
in Fig. 5. 
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Figure 7: Kinematic distributions in events produced in the GDM benchmark scenario r], as 
in Fig. 5. 
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• Gravitino: We consider the gravitino to be observable in GDM scenarios where the 
NLSP can be stopped and its decays observed. This was obviously not the case in scenario 
5, considered above, and perhaps not in scenarios rj either, as we discuss later. On the 
other hand (see below), it seems possible to obtain a substantial sample of fi decays in 
scenario e, so we consider the gravitino to be indirectly observable in this case. 

• Sleptons: The fi has a distinctive time-of-flight signature, and we consider that it 
could be detected with some efficiency in both ATLAS and CMS in all the scenarios e, C, ?7- 
As seen in panels (d) of Figs. 5, 6 and 7, it is generally produced quite centrally and in 
association with a considerable number of high-£^r jets and/or leptons. As discussed below, 
its mass can probably be measured with an accuracy ~ 1 % in all three scenarios, and a 
significant sample of decays of stopped fi's should enable its lifetime to be measured at 
point e. Even in scenarios expects a sample of several hundred events with qjiQu 
production followed by qr — > q{x — tit) decay on one side, and qr — > q{x — ^r^) decay 
on the other side, followed by Ir — > ifiT decay. We therefore expect the £r to also be 
observable at all three points. In scenario e, one should also be able to reconstruct the 
cascade Ql ^ ?X2 followed by X2 ^ ^l^ and then £l ^ X^- Knowing already the x mass 
from the analysis of qR decays, one should be able to reconstruct the II and xi masses at 
point e. At points ?7, we expect a sample of a few dozen events, which might be sufficient 
to argue that > > > nif^, and thereby provide some discrimination against 
GMSB models [43] but would be insufficient to reconstruct the heavier slepton masses. 
In summary, therefore, we consider all the charged sleptons to be observable at point e, but 
only fi and £r at points (^,7]. 

• Sneutrinos: We do not consider these to be observable in any of the the three scenarios. 

• Gauginos: The lightest neutralino x should be observable in all three scenarios e, C, ?7, 
as a resonance in fiT combinations, for example in qr — > q{x t^it) cascade decays which 
have branching ratios of 92/75/69 % in the three scenarios. On the other hand, the second 
neutralino X2 is probably observable only in ql decays, and only in scenario e. We do not 
consider the charginos and heavier neutralinos to be observable at any of these points. 

• Higgs bosons: The h should be observable in the three scenarios e to 77, but the heavier 
Higgs bosons are not expected to be observable in any of them. 

• Squarks: At point e, all squark flavours, excluding the {2 but including the ii and 

^1,2 (which appear in 10/15/11 % of ^ decays), should be observable. The spartners of the 

u, d, s, c quarks could also be observed at points C, ?7, but not the ii^2 and 61^2- 

we discuss later, the potential discriminators between gravity-mediated GDM models and GMSB 
models include the sequence of sparticle masses and the pattern of ino mixing. 
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• Gluinos: According to our standard criteria, these should be observed in all three 
scenarios eX,V- 

We have made a first examination of issues in the reconstruction of sparticle cascade 
decays in GDM benchmark scenarios e, r), using as initial building-blocks the final-state 
fi and T. The px distributions for hadronic r-decay jets at the three points are shown in 
Fig. 8(a, c, e), where we see that those at points (^,r) are significantly harder than at point 
e. The mis-tag probabilities we take from [26] are adequate for identifying the large-p^ t 
hadronic jets in GDM sparticle-pair-production events, which would have been collected by 
the normal large-pT trigger at the LHC. Using these tagging estimates, and combining the 
candidate r hadronic jets with the fi tracks, which are assumed to be measured with an 
accuracy 5p/p — 5 %, we obtain the r — fi invariant-mass distributions shown in panels (b, 
d, f) of Fig. 8, respectively. In each case, we see a clear signal due to x — > rfi decays. 

We have then considered the reconstruction of higher cascade decays in scenarios e and 
the latter being very similar to point rj. We find peaks in fi—T—i combinations corresponding 
to the £l,r, and in fi — t — i — i combinations corresponding to the X2- However, the 
combinatorial backgrouds have shapes quite similar to the signals. Full studies of these 
scenarios lie beyond the scope of this survey, but it does seem that sparticle cascades can 
be reconstructed in these scenarios, analogously to what was shown previously for scenarios 
with a X LSP. 

5.2 Detectability at e'^e~ Linear Colliders 

As in our previous studies [7, 9] , our criteria for the observability of supersymmetric particles 
at linear colliders are based on their pair-production cross sections. 

• Particles with cross sections in excess of 0.1 fb are considered as observable, because 
they would give rise to more than 100 events with an integrated luminosity of 1 ab~^. 

• The lightest neutralino x is considered to be observable only through the decays of 
heavier supersymmetric particles. 

• Sneutrinos are considered to be detectable when the sum of the branching fractions for 
decays which lead to clean experimental signatures, such as i>i —>■ x^£^ {i = e, /i, r) and 
Uj. — > W~^fi~, exceeds 15 %. 

• The 77 collider option at a linear collider would allow single production of heavy neutral 
Higgs bosons via the s-channel processes 77 — > A and 77 — > H, extending the reach up to 
375 GeV for 0.25 TeV e± beams, 750 GeV for 0.5 TeV e± beams, 2.0 TeV for 1.5 TeV e± 
beams and 3.75 TeV for 2.5 TeV beams. A 77 collider may also be used to look for 
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Figure 8: The t pt distributions and the r — fi invariant-mass distributions for (a, b) 
benchmark scenario e, (c, d) benchmark scenario ( and (e, f) benchmark scenario rj. 
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gluinos [42], but we do not include this possibility in our analysis. 

• Finally, we assume that a metastable fi could be detected at any linear e+e~ coUider 
with more than 100 events, and note that the mass could be measured more accurately than 
at the LHC, by measuring the production threshold as well as (see the next Section for 
further discussion). 

As previously [7, 9], we consider e+e~ collision energies = 0.5 TeV, 1.0 TeV, 3 TeV 
and 5 TeV, and also the combined capabilities of the LHC and a 1-TeV hnear collider. 

• For ^/s — 0.5 TeV, at NUHM benchmark point a the x, X2: xt ^^e lighter sleptons 
i^Rj^R, T^i would be observable. The prospects at benchmark point (3 are similar, except that 
the would not be observable. At benchmark point 7, all the inos except the xf would 
be observable, and the sleptons jiR, cr, f\ would be observable in X4 decays. The lightest 
h is always observable and the if, A can be produced in 77 collisions for these benchmark 
points. 

The prospects for the GDM benchmark points are not so good: apart from /i, only 
li-Ri ^L,R, 7^1,2 are observable and then only in the low-mass scenario e that was chosen at 
the tip of the low-mo GDM wedge in mSUGRA parameter space. No squarks or gluinos are 
observable in any NUHM or GDM scenario at = 0.5 TeV. 

• For a/s — 1 TeV, all the neutralinos, charginos and sleptons (both charged and neutral) 
become observable in scenarios a, (5 and 7, and the heavy Higgs bosons can now be pair 
produced in e'^e" collisions directly. In the GDM scenario 5, associated XX2 production 
becomes observable, albeit with a low event rate, and jiR, cr, fi production can be detected 
in X2 decays. At point e, associated XX2 production should again be observable, and probably 
also associated XXSA production, as well as x^ pa-ir production, associated x^X2 production 
and pair-production of all the charged sleptons. Finally, at points C and 77, only x, jiR, cr 
and fi are expected to be observable. 

• For ^/s — 3 TeV, all the Higgs bosons, neutrahnos, charginos, sleptons and squarks 
would be observable in each of the scenarios a, /3, 7, e. The same is true at benchmark 
point 5, with the exception of the left-handed first- and second-generation squarks, and the 
right-handed squarks ^2 and 62, because of their low rates. At benchmark points ^,77, one 
should observe all the weakly-interacting sparticles, but (with the exception of the ii) the 
squarks would still be out of kinematic reach. At benchmark point 7, also the gluino will be 
observable in squark decays. 

• For a/s — 5 TeV, the story is simple: all the sparticles except g would be observable in 
all the proposed NUHM and GDM benchmarks, and also the g at point 7. 

^^The heavier sleptons are visible via associated enCL and fif2 production. 
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5.3 Summary 

Fig. 9 summarizes the numbers of different species of MSSM particles visible at different 
accelerators. We see that the LHC provides good coverage for strongly-interacting sparticles 
(uppermost green and pink bars) in all the NUHM and GDM benchmarks scenarios con- 
sidered, whereas its coverage for weakly-interacting sparticles (middle red and blue bars) is 
rather uneven. The lightest Higgs boson is always detectable (lowest light blue bars) and 
in some cases also heavier Higgs bosons. In scenarios 7,e, a linear e^e~ collider with 
y/s = 0.5 TeV would provide useful extra information about some weakly-interacting spar- 
ticles. In all cases, it would provide detailed measurements of one or more Higgs bosons. A 
linear e+e^ collider with ^/s = 1.0 TeV would provide better information on both weakly- 
interacting sparticles and Higgs bosons, but still no information on squarks or gluinos. The 
combination of the LHC and a 1.0-TeV linear collider would provide good coverage overall, 
but this would still be incomplete in scenarios S,(,Vj particular. A linear e~^e~ collider 
such as CLIC with = 3.0 TeV would provide detailed studies of all the weakly-interacting 
sparticles and Higgs bosons in all the scenarios studied, and also provide new opportunities 
to study squarks in scenarios a, j3,'~f,6,e. Finally, CLIC at 5.0 TeV would provide detailed 
measurements of all the MSSM particles except possibly the gluino, for which one would still 
rely on the LHC 

6 The Stau NLSP in GDM Scenarios 
6.1 Production and Detectability at the LHC 

In the mSUGRA GDM models studied here, all supersymmetric events yield a pair of fi 
NLSPs. The astrophysical BBN/CMB constraint prevents the fi lifetime from exceeding 
~ 3 X 10^ s [34], and we do not discuss here fi lifetimes smaller than 10'^ s. Charged NLSPs 
with lifetimes in this range would appear to a generic collider detector like massive stable 
charged particles, and the three benchmark scenarios e, rj studied here span this range of 
lifetimes. 

Fig. 10 shows the distributions of the non-relativistic factor f3'~f expected for the fi's in 
these mSUGRA GDM scenarios from cascade decays of squarks and gluinos at the LHC. The 
great majority of the fi's produced at the LHC are far from being ultra- relativistic, and so 
should yield exotic time-of- flight (TOF) and/or dE/dx signals The same would be true 

-"^^Unless one could also observe 77 gg [42], a possibility not considered here. 

^°We are not optimistic about the prospects of detecting these signals in hadron-coUider events without 
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Figure 9: Summary of the numbers and types of MSSM particles that may be detected at 
various accelerators in the NUHM and GDM benchmark scenarios considered here. The 
gravitino may also be indirectly 'observable' at point e. We emphasize that these numbers 
are estimates that need to be validated by experimental simulations. As in [7, 9], we see that 
the capabilities of the LHC and of linear e^e~ colliders are largely complementary. We re- 
emphasize that mass and coupling measurements at e^e~ colliders are usually much cleaner 
and more precise than at the LHC, where, for example, it is not known how to distinguish 
the light squark flavours. 
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of fi's produced at the ILC, if its centre-of-mass energy reaches above the pair-production 
threshold in the corresponding scenarios, namely 310, 680 and 650 GeV in benchmarks e, ( 
and 77, respectively. CLIC with a centre-of-mass energy of 3 TeV would be able to produce 
fi pairs with masses < 1.5 TeV, and hence probe the GDM wedge out to mi/2 ~ 4.5 TeV. 

The key signal for GDM with a fi NLSP would in general be the coincident appearance 
in adjacent or nearby bunch crossings of generic high-pT triggers and subsequent 'muon' 
triggers. Such coincidences would be very rare for conventional trigger rates, which would 
give coincidence rates < 10"^ in adjacent bunch crossings. 

The fi NLSP is often sufficiently non-relativistic, /37 < 1, that it would not exit an LHC 
detector such as ATLAS or CMS before the next bunch crossing (25 ns after the event in 
which the fi was produced), in which case its tracking information might be lost. This 
issue could be addressed by reading out of the detector all the tracking signals that occurred 
within several crossing times following an 'interesting' event. 'Interest' would normally be 
defined by a conventional high-pr lepton or calorimetric trigger. As seen in Figs. 5, 6 and 
7, most supersymmetric events would indeed pass the normal ATLAS and CMS criteria for 
'interesting' events. 

If a sample of interesting candidate events can be identified, one possible fi search strategy 
would be to select out of the usual high-p^ lepton and calorimetric triggers a subsample of 
events suspected of containing fi NLSPs. Even if the muon systems do not trigger on the 
fi's, the muon drift tubes of both ATLAS and CMS integrate signals over a number of bunch 
crossings, such that hits of particles which are out of the normal 25-ns time window can still 
be recorded with the event. One would, however, need to adapt the track reconstruction 
software so as to allow for such signal time shifts and recuperate the full information. At the 
moment, in the absence of any good reason to take seriously such scenarios with massive, 
slow-moving charged particles, the only experimental strategy required from the ATLAS and 
CMS Collaborations is to avoid precluding the possibility of such a buffered readout, should 
it subsequently appear worthwhile for searches in GDM or other scenarios. 

Alternatively, one could use the presence of a high-p^ charged particle in the muon 
system as a primary trigger, and then look back through earlier bunch crossings for evidence 
of other high-p^ jets and/or leptons, that would already have triggered the detector and 
been recorded. 

other distinguishing features, such as the Drell-Yan production of n pairs. 
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6.2 Measuring the Stau Mass 

A crude estimate of the obtainable fi mass resolution can be derived by propagating the un- 
certainties in the momentum measurement Ap and in the TOF resolution At, as determined 
in a detector at a distance L: 

For ATLAS and CMS, the expectations are Ap/p ~ 1 — 10 % and At ~ 1 ns at a distance 
of ~ 5 m. Since the peak value of /37^, is ~ 2, as seen in Fig. 10, we estimate that in each 
event 

^ = (0.01-0.10)00.12. (7) 

We therefore estimate that the fi mass could be measured with an error of 10 - 20 % in 

each event, which could be reduced by selecting low- momentum events as shown in Fig. 11, 

and further by combining measurements in many events. Panel (a) shows the distribution 

of the second term in (6) for fi's produced in a sample of events from benchmark point e, 

and panel (b) shows the distribution of the corresponding mass resolution AM obtainable 

event-by-event. Selecting now a sample of the 10 % of fi's with the lowest values of /37^, and 

hence, according to (6), those with the smallest values of AM and assuming that Ap/p 

has a Gaussian error of 5 %, we obtain panel (c) of Fig. 11. We see that almost all these 

individual events have AM/M < 10 %. The same is true for an anologous sample of events 

produced in a simulation of benchmark point C, as seen in panel (d) of Fig. 11 and also for 

point 77 (not shown). Therefore, if one could obtain a total sample of ~ 1000 fi events, and 

if systematic effects and correlations could be controlled, combining the best-measured 100 

events could yield AM/M < 1 %. 

To help assess the importance of such a measurement, we display in Fig. 12 how a 1 % 

fi mass measurement could be combined with a determination of the supersymmetric mass 

scale mi/2 to determine the allowed range of parameters in the (mi/2,mo) plane. If mi/2 

were known exactly, the error in ruf^ would correspond to 5mQ ~ 20 GeV, which should 

be compounded with the error induced by propagating the uncertainty in mi/2. Such a 

^^In contrast, in the context of gauge-mediated supersymmetry-breaking models (GMSB) [43] with masses 
similar to those in the GDM models considered here, Ref. [44] considered the measurement of particles 
in the upper range of (3 before the next bunch crossing. This would not be optimal for measuring the 
metastable particle mass, but it was nevertheless estimated that a precision AM/M < 1 % could be attained. 
Measuring the lifetime of the metastable particle inside the collider detector was also considered in such 
GMSB models [44], but in GDM models this would be feasible only after first stopping the NLSP, since its 
lifetime is much longer, as we discuss below. 
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Figure 11: (a) The TOP contribution to the mass resolution AM/M (6) in the stau NLSP 
mass measurement, (b) the full uncertainty AM/M found for the complete event sample of 
benchmark point e, (c) AM/M for the 10 % of events at point e with the lowest values of 
(3'-)'^, and (d) the same distribution for point ( (the result for point rj is very similar). 
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determination of mo would enable a ballpark estimate of the fi NLSP lifetime to be made, 
within this mSUGRA GDM framework, enabhng a strategy to search for its decays to be 
better focused, as also seen in Fig. 12. 

The value of mi/2 could perhaps be determined by measuring the gluino mass, but here 
we discuss the use of the total supersymmetric cross section. Comparing the total sparticle 
production cross section in scenario e with those for scenarios ^,77, we see that atot ~ "^1/2' 
approximately. The statistical error in measuring atot for either of the scenarios ^,77 would 
be about 1.5 %, but we expect that the systematic and theoretical errors would be larger. 
Neglecting theoretical errors, an experimental error of 5 % would enable mi/2 to be estimated 
with an uncertainty < 1 % from the total cross section alone: see the vertical shaded band in 
Fig. 12 Assuming also a measurement uncertainty of ~ 1 % in m^^, we see from Table 2 
and Fig. 12 (diagonal shaded bands) that this should be sufficient to distinguish between 
scenarios V at the 5-a level. However, this discrimination would be lost if the error in 
either mf^ or mi/2 rose to ~ 5 %. 

One could in principle also look for slowly-moving massive charged particles via an anoma- 
lous dE/dx signal, for instance using high-threshold signals in the ATLAS TRT detector, 
and maybe using the electromagnetic calorimeter of either experiment A good dE/dx 
measurement is a design feature that could be considered for the central trackers of future 
detectors at the ILC and/or CLIC. 

6.3 Distinguishing GDM Benchmarks from GMSB Models 

We now consider the possibility of using spectroscopic measurements at the LHC to distin- 
guish the GDM benchmark scenarios considered here from minimal gauge-mediated models 
of supersymmetry breaking (mGMSB) We first consider the spectroscopic properties of 
the 'easy' GDM scenario e, and compare them with mGMSB models with the same value of 
mi/2. The principal discriminants we consider are the masses of the %, X2, the fact that the 
£r is significantly lighter than the x, and the qn mass. 

We recall that mGMSB mass spectra are characterized typically by the messenger index 
N, a mass A that sets the overall gaugino mass scale: Ma = {aa/A'n)N K (a = 1, 2, 3) and the 

^^A more conservative error estimate of ^ 25 % woiild yield the wider band indicated by dashed vertical 
lines. 

^•^It might also be interesting to add to ATLAS or CMS a specialized time-of-flight detector, either in the 
cavern itself or suspended in the access pit outside it [45] . 

^^LHC measurements of supersymmetric cascade decay branching ratios might also help discriminate, but 
we do not consider them here. Since the NLSP lifetimes are very different at our GDM benchmark points 
from mGMSB models, the detection of decays into gravitinos, discussed in the next subsection, would also 
help in the discrimination. The ILC would be able to distinguish GDM from mGMSB very easily. 
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Figure 12: The potential impacts of prospective measurement errors of 1 % and 5 % for 
nif^ (diagonal hands and lines) and mi/2 (vertical hands and lines) as constraints in the 
(mi/2,'n^o) plane for GDM models in the mSUGRA framework. The smaller errors would 
enahle the henchmark scenarios ( and rj to he distinguished, and the possihle NLSP lifetime 
to he estimated. The thin solid lines are laheled hy the logarithm of the NLSP lifetime in 
seconds. 
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Figure 13: Sparticle masses in minimal GMSB models with N = 2 and (a) mi/2 = 440 GeV, 
tan/3 = 15 and (b) mi/2 = 1000 GeV and tan/5 = 22, as functions of A for (a) 
M = 800,300,200 TeV (diagonal dashed, solid and dashed lines, respectively), and (h) 
M = 8000,1000,300 TeV (diagonal dashed, solid and dashed lines, respectively), compared 
with GDM predictions (horizontal dot- dashed lines). The best agreement between the two 
models is found for (^aj A ~ 67 TeV and (b) A c:^ 155 TeV, respectively. 

scalar masses: mg. oc A^A^, and an input scale M from which the RGEs are used to evolve 
the sparticle masses down to Qewsb = 2 TeV. 

For A^ = 1, the in is always heavier than the x, for A^ = 2 the x ~ mass difference 
may be approximately the same as at the GDM benchmarks if M ~ 200 — 300 TeV, and 
for A^ = 3 the mass difference is significantly larger. Therefore, we can discard A^ = 1 and 
concentrate on = 2 while retaining = 3 as a second option. As seen in Fig. 13(a), a 
GMSB model with N = 2 has approximately the same values of and m^^ the GDM 
benchmarks if A ~ 67 TeV, whereas the best value of A would be somewhat lower for = 3. 
In each case, rrig is very similar in the GDM and mGMSB model for the best value of A. 
However, in the N = 2 case rriq^ > rrig, whereas in the = 3 case rriq^ < rrig. Thus, an 
LHC measurement of — m^^ has the potential to exclude mGMSB with A^ = 3 and that 
of rrig — rriq^ to exclude N = 2. 

We have also considered the 'average' spectroscopic properties of GDM scenarios ( and 
rj, and compared them with mGMSB models with the same value of mi/2 = 1000 GeV and 



39 



tan (3 — 22. The principal discriminants we consider are again the masses of the x, X2, (-r and 
qr. As in the case mi/2 = 440 GeV, we again find that N — 1 gives a iu heavier than the x, 
and hence can be discarded, whereas N — 3 gives too large a mass difference — fni^. On 
the other hand, as seen in panel (b) of Fig. 1?,, N — 2 gives a qr heavier than the g whereas 
N — 3 again gives a lighter qr. As in the previous case, an LHC measurement of — 
has the potential to exclude mGMSB with N — 3 and that of rrig — niqj^ to exclude N — 2. 

We conclude that a combination of these mass measurements and the different lifetimes 
of NLSP decays would be sufficient to distinguish the GDM models of the type considered 
here from mGMSB models. 

6.4 Stau Trapping and the Detection of Decays 

If produced fi's arc sufficiently slow-moving, they may be stopped inside the detector or 
its neighbourhood. We consider three possibilities: that the fi may be trapped inside the 
detector itself, or in adjacent water tank or calorimetric detector, or in the walls of the 
experimental cavern. In the case of the LHC, the trapping rate can be calculated using the 
spectra shown in Fig. 10 and the known rates of energy loss by charged particles passing 
though different types of matter [46]. As representative examples, we consider Iron - as in 
an experimental calorimeter - and Carbon - which has similar stopping power to water and 
less than any other plausible detector and/or surrounding material. 

We display in Table 7 the numbers of fi's expected to be produced in GDM benchmark 
scenarios eX,V at the LHC with 100 fb~^ of integrated luminosity, using the cross sections 
shown in Table 6 calculated with PROSPINO [40]. We consider numbers of particles produced 
with /37 < 0.25(0.5) and show in each case the corresponding ranges in Carbon and Iron. 



Model 


e 


c 


V 


Number of particles with 


850 


7 


7 


/57 < 0.25 








Range in C (cm) 


60 


136 


129 


Range in Fe (cm) 


29 


65 


61 


Number of particles with 


7700 


100 


90 


/?7 < 0.5 








Range in C (cm) 


600 


1360 


1290 


Range in Fe (cm) 


290 


650 


610 



Table 7: Numbers of slow-moving ti 's produced with 100/fb at the LHC in GDM benchmark 
scenarios eX,V) the corresponding ranges in Carbon and Iron. 

Benchmark scenario e is in a different class from scenarios (,r). The numbers of slow- 
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moving particles are orders of magnitude larger, and secondly the ranges are shorter typically 
by a factor of two. Both these features are directly related to the sparticle mass scale, which 
is set essentially by mi/2. In the case of benchmark e, hundreds of fi's would be trapped 
within either the ATLAS or CMS calorimeter, and thousands more would be trapped within 
a few metres of surrounding material. On the other hand, in the cases of scenarios C, ?7, only 
a handful of particles would be trapped within either detector, and only a few dozen events 
would be trapped within ~ 10 m of surrounding material. 

The situation at the ILC would be very favourable for any of the scenarios considered. 
The above reference values of (3^ correspond to ^/s — 2mf^ x 1.0625(1.25), and there should 
be no problem tuning the beam energy very close to the fi mass and obtaining large samples 
of fi's stopped within the calorimeter. The same would be true at CLIC for mSUGRA GDM 
models with mi/2 < 4.5 TeV. 

In the case of the LHC, unfortunately there is very little room left in the ATLAS cavern 
for a trapping water tank or calorimeter, and it would be difficult to envisage inserting any 
detector over a metre thick between the barrel and the cavern walls. On the other hand, in the 
case of CMS there may be some more room after restructuring the infrastructure (balconies 
and services), permitting the installation of an O(kton) trapping detector, as discussed 
in [47]. There would be more room in the forward direction at CMS, but this possibihty 
would have limited angular acceptance. Moreover, the fi pseudorapidity distributions are 
generally very central, within the ATLAS or CMS acceptance, as shown in Figs. 5, 6, 7, so a 
forward trap would not be very efficient. In the case of the ILC or CLIC, if the experiments 
at the LHC reveal interest, it would be possible to design the experimental areas ahead of 
time so as to allow for a trapping detector. 

In the interim, we speculate on the alternative possibility of looking for the decays of 
fi's that are trapped in the walls of the ATLAS and CMS caverns. One possible strategy 
would be to use the tracking information from CMS or ATLAS to determine the fi's impact 
point and angle, then bore a hole into the wall and extract a core with an optimal chance 
of containing a trapped fi. The tracking systems of CMS and ATLAS should each yield an 
experimental uncertainty in the impact point that is about half a cm, and a corresponding 
angular error ~ 10~^ radians. Using the standard formula [46] 



for the 98 % C.L. width of the projected distribution of the multiple scattering angle, where 
P,p are the fi velocity and momentum and x the penetration depth relative to the scattering 
length Xq, we find typical values < 10~^, within the experimental angular error. As we 
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can see from Table 7, one might want to extract O(IOO) to 0(10000) 'interesting' cores with 
dimensions ~ 1 cm xl cm xlO m each year. This technique might be appropriate for the 
upper part of the mSUGRA wedge shown in Fig. 1, where the fi hfetime is measured in 
weeks, such as scenarios e, However, this is unhkely to be feasible in the lower part of 
the mSUGRA wedge, e.g., at point 77, because radiation levels in the LHC caverns would 
preclude access on the necessary time-scale Tf^ ~ 5 hours. 

The baseline operating plan for the LHC foresees one multi-month stop each winter, and 
half-a-dozen two-day technical stops at regular intervals during the rest of the year. Each 
of these would provide an opportunity to extract a limited number of cores from the cavern 
walls. This would be interesting if the fi lifetime is several weeks or more, as in benchmarks 
e, C, but not point 77. 

We have also considered the possibility of measuring directly the mass of a stopped fi in 
a mass spectrometer. A typical extracted core of size 1 cm x 1 cm x 10 m would contain 
~ 1 X 10^^ protons. On the basis of estimates of the mass of the fi and the velocity of the 
specific particle being sought in the core, we estimate that a 'high- interest' sample of about 
10 % of the length of each core might be selected for exploration in more detail using a mass 
spectrometer. For comparison, we note that the best available upper limit on the relative 
abundance in water of a positively-charged stable relic particle with mass between 40 and 
400 GeV is 1 X 10~^^ [48]. It might therefore be feasible to pass the high-interest samples of 
each core through a mass spectrometer and measure the fi mass very precisely. The issue 
would be how quickly this study could be completed, in comparison with the fi lifetime, as 
discussed above. 

Another possibility might be to look for upward- or sideways-going muons coming out of 
the wall, produced by fi decays in the neighbouring rock. We estimate that typical muon 
momenta would be tens of GeV, in which case they should be able to traverse tens of metres 
of rock. However, the acceptance for decays back into the cavern would not be large unless 
the fi decays within a few metres of the cavern wall. In the benchmarks studied, detecting 
these might be feasible for the thousands of fi's produced with < 0.5 in scenario e, but 
looks very marginal for the few dozen fi's produced with < 0.5 in scenarios ^,77, which 
would also have longer ranges, diminishing the angular acceptance for the decays. 

We have compared the possible measurement of the 'albedo' due to fi — > r — > // decay, 
which has a branching ratio of ~ 16 %, with the irreducible background due to the known 
atmospheric u ^ /i flux. Assuming that the gravitino mass is small compared with nif^, 
the characteristic /i energy will be ~ nif^/G, corresponding to ~ 25 GeV for point e and 
~ 50 GeV for points e, (. If we consider a representative LHC or ILC detector with hnear 
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dimensions 20 m x 20 m and fi decay at a characteristic distance ~ 10 m, the detector 
subtends 1/6 of the total sohd angle, namely ^ steradians. However, fi decays in the upper 
hemisphere surrounding the detector will surely be drowned in cosmic-ray /i background, 
so we consider only fi decays in the lower hemisphere. Therefore, only 1/12 of the decay 
muons are in principle observable, corresponding to just 1.3 % of the stopped fi decays. As 
seen from Table 7, at the LHC this would give O(IOO) events in benchmark scenario e, but 
only 0(1) event in either of scenarios rj. For comparison, the MACRO experiment [49] has 
reported a sample of about 900 upward-going atmospheric u ^ /i events passing through a 
detector of area 76 m x 12 m in five years of operation, corresponding to ~ 80 events/year 
through our nominal 20 m x 20 m collider detector. The energy spectrum of the MACRO 
'up-through' II sample has a broad peak around 50 GeV, so there is no clear separation in 
energy between the fi — > r — > // signal and the atmospheric — > // background. On the 
other hand, the tracking system defines the direction in which a candidate fi exited the 
detector, and the momentum measurement constrains the distance at which it is likely to 
have stopped. Together, these measurements define the direction from which a candidate 
'albedo' decay muon might emerge from the wall into any part of the detector. Using this 
information, it might be possible to detect the 'albedo' decay muons at the LHC in scenario 
e, though not in scenarios V- the other hand, detection at the ILC could be optimized 
with a dedicated detector. 

Comparing the mSUGRA spectra for benchmark scenarios e, C, 77 with the sample of 
GMSB models in [43], we see that there exist GMSB models with identical values of (mf^, {rriq)) 
We are therefore pessimistic that these spectroscopic measurements at the LHC will be 
enough alone to distinguish mSUGRA from GMSB. 

We have also considered the possibility of using fi — > tG decay kinematics to constrain 
directly the mass of the gravitino G, and hence perhaps also distinguish between mSUGRA 
and GMSB scenarios. In the three fi NLSP scenarios e, C, V considered here, the mean r 
energies are (Er) — 74, 155, 160 GeV, respectively. On the other hand, in GMSB scenarios 
with the same values of mf^, one would have (E-,-) — 75,170,161 GeV, respectively. We 
recall that rrif^ would be measurable with an accuracy that is probably not much better 
than 1 %, which already removes any chance of measuring the percentage difference in (Et) 
in scenarios e, 77, where it is 0{1) %. On the other hand, the percentage difference in scenario 
C is about 10 %, which should be measurable in principle with enough events and accurate 
energy measurements. However, in practice, as we have discussed above, a sample size of 
more than 100 fi decays, as would be required for a 10 % measurement of (Ej-) in scenario 
even if Er could be measured perfectly event-by-event, could not be obtained inside ATLAS 
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or CMS with the 100 fb~^ integrated luminosity assumed here. 

On the other hand, at the ILC nif^ should be measurable with an accuracy better than 
0.1 %, and it should be possible to obtain a large enough sample to measure (Er) with an 
accuracy of a few %. We therefore think that distinguishing scenario ( from a GMSB model 
with the same nif^ should be possible, but making the same distinction for scenarios V 
would be much more challenging, in view of the very accurate (Ej-) measurement required. 

7 Conclusions 

We have discussed in this paper various alternatives to the constrained MSSM (CMSSM) 
framework with universal soft supersymmetry breaking, which has been used in most pre- 
vious benchmark studies of supersymmetric signatures at the LHC, ILC and CLIC. Specifi- 
cally, we have considered both a less constrained framework with non-universal Higgs masses 
(NUHM) and a more constrained gravitino dark matter (GDM) framework inspired by min- 
imal supergravity (mSUGRA). As we have shown, both of these scenarios offer distinctive 
phenomenological signatures that were absent in the previous CMSSM studies, at least when 
cosmological constraints from WMAP and other experiments were taken into account. For 
example, Gravitino Dark Matter (GDM) becomes a generic possibility. 

In the case of NUHM models, the freedom to vary the scalar masses relative to the gaugino 
mass mi/2 opens up the possibility of different characteristic heavier neutralino decays such 
8^ X2 — ^ Xh, X-^, which were disfavoured along the WMAP hues allowed in the CMSSM. 
Three of the new benchmark scenarios we have proposed and examined here are of this type. 
One has a dominant X2 — decay mode (point a), one has a dominant X2 — ^ X^ decay 
mode (point (3), and one has non-resonant X2 X^^^" decays (point 7). The h signal at 
point /3 should easily be detected at the LHC, and a fortiori at the ILC or CLIC. These 
features open up new possibilities for reconstructing the cascade decays of heavier sparticles 
and measuring better the mass of the lightest neutralino x- 

In the case of GDM models, there are two generic possibilities for the next-to-hghtest 
supersymmetric particle (NLSP). Either it is neutral and presumably the lightest neutralino 
X, as exemplified by benchmark scenario 5, or it is charged and presumably the lighter stau 
slepton fi, as exemplified by benchmark scenarios e, C and 77. 

In the former case, which is a CMSSM scenario with the gravitino specifically chosen 
to be fight, the colfider signatures are very similar to those for scenarios with a stable x, 
since generic decays of the unstable x occur far outside the collider detector. However, 
the sparticle spectrum is somewhat different from what it would be along a WMAP line. 
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providing an opportunity to distinguish this possibihty, for example in ^ — > 6& decays. 

The three benchmark scenarios e, ( and 77 with a fi NLSP span the range of possibihties 
that are hkely to be detectable at the LHC or ILC within the mSUGRA framework. Point 
e was chosen to have sparticles as light as possible, with heavier spectra at points C, 77. At 
point C the fi lifetime is close to its theoretical maximum, and at point 77 the fi lifetime is 
close to the minimum we consider. 

In all three cases, the fi should be detectable as a slowly-moving metastable particle, and 
TOP measurements at the LHC might enable its mass to be measured with an accuracy at 
the % level. We have discussed how these mSUGRA GDM scenarios could be distinguished 
from mGMSB models. 

As pointed out in [35, 36], a small fraction of the produced fi's may be stopped inside 
the collider detector and/or in surrounding material, either cavern walls or a specialized 
stopping detector. At the LHC, there is hmited space in the constructed caverns for such 
a specialized detector [47] , but this could be envisaged in designing intersection regions for 
the ILC or CLIC, if needed. We have discussed a couple of strategies for detecting fi's that 
might have stopped in the cavern walls surrounding ATLAS or CMS, either by extracting 
cores from the surrounding wall material or by looking for muons from the decays fi — > t+ 
gravitino followed by r — > // + fgP^. The validities of these strategies depend on practical 
feasibility as well as the fi lifetime, which is uncertain by two orders of magnitude even in 
the restricted mSUGRA framework discussed here. 

Our preliminary investigations indicate that the LHC has good capabilities to discover 
supersymmetry in each of the benchmark scenarios proposed, and may be able to uncover 
several different sparticles in each case. However, in general the LHC experiments would 
not be able to make many very accurate spectroscopic measurements. In several of the 
benchmark scenarios, an ILC with a/s — 0.5 TeV would already be able to make several 
precise measurements of weakly-interacting sparticles, and a 1 TeV ILC would extend this 
capability to all the benchmarks studied here. In all the scenarios studied, the additional 
measurements possible with CLIC at ^/s = 3 or 5 TeV would contribute significant added 
value to the exploration of the new physics previously uncovered by the LHC and followed 
up by the ILC. 

These new benchmark scenarios merely scratch the surface of the myriad possibilities 
open in supersymmetric phenomenology once one explores the high-dimensional parameter 
space of soft supersymmetry breaking. These examples should offer general encouragement to 
think outside the CMSSM box, and suggest the likehhood that even more exotic possibihties 
might be waiting out there, perhaps including one chosen by Nature. 
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